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The possibility of creating a Bose condensate of trapped atoms in a non-ground state is suggested. 
Such a nonequilibrium Bose condensate can be formed if one, first, obtains the conventional Bose 
condensate in the ground state and then transfers the condensed atoms to a non-ground state by 
\f^ ' means of a resonance pumping. The properties of ground and non-ground states are compared and 

plausible applications of such nonequilibrium condensates are discussed. 
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I. INTRODUCTION 



The recent realization of Bose-Einstcin condensation of dilute atomic gases in magnetic traps [1-3] has opened 
a rapidly expanding field of studies of condensate properties. There has been a splash of both experimental and 
theoretical activity on this subject (see books [4,5]). 

Atoms trapped in a confining potential possess the discrete spectrum of states. At high temperatures these states 

fj • are occupied, according to the quantum Bose-Einstein distribution, so that no state is occupied macroscopically. 

I— ' ' Under the macroscopic occupation of a state one implies that the number of atoms in this state is proportional to the 

total number of atoms in the system. An important consequence of quantum statistics is that, when the system is 

K^ ' cooled down below some critical temperature, bosons pile up in the lowest energy state of a confining potential. The 

Qv 1 macroscopic population of the quantum-mechanical ground state of a confining potential is the characteristic feature 

\^ ' of Bose-Einstein condensation. 

f^ A natural question that can be raised is: Is it possible to realize the macroscopic population of some other quantum- 

OQ mechanical state rather than the ground state, or in addition to the latter? That is, can one produce a Bose condensate 

in a non-ground state? The answer to this question is interesting by its own. And, if that is possible, several new 
important applications can be suggested. 

For example, recently two overlapping *^i?6 condensates in two different ground state hyperfine levels were created 
"ti ■ [6]. However, for other atoms the simultaneous creation of two condensates in different ground-state hyperfine levels 
r^ ' may be difficult or not feasible [7]. Then the alternative could be the creation of two condensates, one in the ground- 
'T' , state and another in a non-ground-state level. 
'"^ ■ Realizing the macroscopic population of a non-ground state of a confining potential could be a way for producing 

^ [ various spatial distributions in the system of coherent atoms. This may find application for atom lasers for which the 
O ' creation of coherent atomic beams with different spatial modes may be required. 

One more possibility of employing such a non-ground-state condensate could be for studying relaxation processes 

in the quantum degenerate regime. The macroscopic occupation of an excited state would result, as is clear, in a 

nonequilibrium sample. When external forces supporting this state are switched off, the system will relax returning 

to the equilibrium by repopulating the discrete levels of the potential. The observation of this process of relaxation 

C^ in the quantum degenerate regime can provide useful information about this form of quantum matter. 

Finally, when a new system with unusual features is explored, there is always the chance of finding something 
completely unexpected. 

In this paper we describe a possible way of transferring the macroscopic number of atoms from the conventional 
ground-state condensate to a non-ground-state level of the confining potential, thus creating a non-ground-state 
condensate. 



II. RESONANCE PUMPING 

Assume that the Bose gas of neutral atoms has been cooled down so that all atoms are in a coherent condensate state. 
The latter is described by the nonlinear Schrodingcr equation which is often called the Gross Ginsburg Pitaevskii 
equation [8-12]. This equation writes 

in^=Hv, (1) 



where the nonlinear Hamiltonian 



contains the nonlinear part 



dt 



H = H{ip) + Vp (2) 



H{^)^-—V^ + U,{r) + A\^\^ (3) 



and, in general, a time-dependent part Vp related to external fields. The atom-atom interaction is modelled by the 
s-wave scattering interaction with the amplitude 

A = {N - l)ATin^ — (4) 

m 

in which N is the number of atoms in the system, a is the s-wave scattering length, and m is the atomic mass. The 
term Uc is a confining potential. The wave function is normalized to unity: {(p,(p) — 1. Assume that at the initial 
time t = 0, all atoms are in the ground state 

(p(7,0) = ^o(r) (5) 

corresponding to the minimal energy level of the eigenvalue problem 

H{(pn)'Pn = Enipn, (6) 

in which n is a multr index enumerating quantum states. The chemical potential is incorporated into the notation of 
the energy levels En- Temperature is assumed to be much lower than the condensation temperature, since only then 
it is possible to condensate almost all atoms in the ground state. Note that the nonlinear Schrodinger equation (1) 
describes coherent states [13]. 

Since the atoms are assumed to be initially condensed in the ground state, to transfer them to higher levels one 
needs to apply an external pumping field that we take in the form 

Vp = V(j) cos Lot. (7) 

Such a field can be realized by a special modulation of the magnetic field producing the trap. As far as our aim is to 
populate a separate energy level, say a particular level p with the energy Ep, we have to choose the frequency of the 
pumping field (7) satisfying some resonance conditions. 
Denote the transition frequencies ujmn by the relation 

'hLOmn = Em — En (8) 

and the detuning from the chosen particular transition frequency uOpQ as 

Ao; = cij — LOpf). (9) 

The first evident resonance condition is that the detuning must be small as compared to the transition frequency Wpo 
corresponding to the transition from the ground state, with the energy Eq, to the chosen particular state with the 
energy Ep, 

^ « 1. (10) 



In addition, it is necessary that the pumping would not influence the neighboring states, that is the detuning must 
satisfy the inequalities 



Ao 



Wp+l,p 



«1, 



Aw 



^p.p— 1 



<1. 



(11) 



The resonance conditions (10) and (11) are necessary but not yet sufficient, if we keep in mind to populate only 
one particular level. This goal can be reached only when the pumping does not force transitions between other states, 
which can be expressed as the inequality 



Auj 



<1 



{m ^ p, n ^ 0). 



(12) 



The conditions (10) and (11) are easy to accomplish making the detuning sufficiently small. The resonance condition 
(12) is more restrictive requiring that the spectrum {£„}, defined by the eigenproblem (6), would not be equidistant, 
as it happens for a simple harmonic oscillator. Really, if that would be the case, then the pumping of atoms from 
the ground state to the chosen particular state would, at the same time, induce transitions from the latter to another 
equidistant state, and from the latter to higher equidistant states, and so on. In such a case, all atoms will be dispersed 
over a number of states making it impossible to get a macroscopic population of any of them. Fortunately, because of 
the nonlinearity, representing atomic interactions, in the Hamiltonian (3), the spectrum {En} is not equidistant even 
when the confining potential Uc is harmonic. In addition, we may include into the confining potential anharmonic 
terms and regulate its spectral characteristics by varying anharmonicity parameters [14-16]. Moreover, as we shall 
show in Sec. Ill, the spectral properties of the nonlinear Hamiltonian (3) may be essentially modified by varying the 
intensity of interactions (4), for which it is sufficient to change the number of atoms N. Therefore, it is always 
possible to prepare the system for which condition (12) holds true. The situation here is similar to the problem of 
inducing resonant electron transitions in an atom. The latter also contains many electron levels, but, because these 
are not equidistant, it is practically always possible to induce a resonant transition between a chosen pair of them [17]. 
The principal difference between the resonant electronic transitions in an atom and atomic transitions in a confining 
potential is that electronic levels are not equidistant because of the hydrogen-type potential, while the interactions 
between electrons do not play essential role. In such a case, the resonant electronic transitions can be treated in a 
linear approximation. As to the atoms in a confining potential, if the latter is harmonic, then the main role of making 
the energy levels non-equidistant is played by the atomic interactions. This makes the problem principally nonlinear 
and forces to deal with complicated nonlinear equations. 

To describe the time evolution of the system, we have to consider the time-dependent nonlinear Schrodinger equation 
(1). Present its solution as an expansion 



'^(^' ^) = X! Cn{t)'fin(r) CXp f -T^nt 



(13) 



in the basis of the stationary states of the eigenproblem (6). Substituting Eq.(13) into Eq.(l), take into account that 
in the double sum 

y^CfcQ(</?myfc,y;y«)exp(iwfc;t) 

k,l 



the main contribution comes from the term 



y^ |cfe|^(93„(/3fe,(/3fe(p„) 



because other terms containing the oscillating factors, being summed up, on average cancel each other. Then from 
Eqs.(l) and (13) we have 



.^dCn 



i?i-^ =^ I VnmCOSUjt+ ^ Ankm\ckf j C„e*' 

fe(#n) 



(14) 



with the matrix elements 



V„-, 



{ipm,V{r)ip„)^ Amkn ^ Aiifirn'Pk.'Pk^n)- 



The solution of Eq.(14) must satisfy the normalization 

n 

Eq.(14) is a set of equations for the functions c„ = c„(t) enumerated by a multi-index n. Let us separate this set 
into the equation for the ground-state function cq, the equation for Cp representing a chosen state with n — p, and 
the equations for all other c^ with k y^ 0,p. Introduce the notion for the population probability 

n, =n,(t)=|c,(t)|2, (15) 

where the index is either j — 0, or j ^ p, or j = k ^ 0,p. Also, define the parameters 

a ^ ^Aopo ^jj |^o(r )n^p(7)|2rf 7 (16) 

and 

(}^^Vop^^l^*,{7)v{7)^p{7)d7 . (17) 

The solution to Eq.(14) can be presented as the sum of a guiding center plus a small oscillating ripple around 
the latter. The equation for the guiding center is obtained from Eq.(14) by averaging its right-hand side over time 
according to the rule ^ /J" F{t)dt with r -^ oo. During this averaging the exponential e*^'^* is treated as a constant 
since it is a slowly varying factor. Actually, at the pure resonance, when Aw -^ 0, this exponential , e*'^'^* -^ 1, is 
exactly one. 

After realizing the described procedures, we obtain from Eq.(14) the system of equations 

dCQ . i r, iAiUt 

— ^ -laUpCa - -j3e Cp, 

— = -tanoCp - 2^ e cq, 

^=0 (MO,p). (18) 

Since the functions Cj are complex, the system of equations in (18) must be completed by another system either for 
the complex conjugate functions c* or for the amplitudes rij — |cjp. The equations for the latter are 



Im {l3e'^'^'4cp) , 



dt ^ 

dn 



^ = /m (r e-*^-*c;co) , (19) 

and duk/dt = when k ^ 0,p. As the initial conditions we have 

co(0) = l, cp(0) = 0, cfe(0) = 0. (20) 

From the last of the equations in (18), together with the initial conditions from (20), it follows that Ck{t) = for 
k ^ 0,p. Therefore, the normalization condition reads 

no{t) + np{t) = 1, (21) 

which demonstrates that the atoms are concentrated in the ground state and p-level, preferentially. 

Eqs.(18) and (19) form a system of complicated nonlinear differential equations. This system could be solved 
by perturbation theory in two limiting cases, when either |a//3| <C 1 or |/3/a| ^ 1. In the intermediate regime, 
when \a/P\ ~ 1, perturbation theory is not applicable. A general solution, valid for arbitrary relation between the 
parameters a and /?, can be obtained by employing the method of scale separation [18-20]. This can be done by 
noticing that the functions Cq and Cp contain time-dependent imaginary factors absent in ng = |cop and Up = |cpp. 



that is, the time variation of cq and Cp is faster than that of no and Up. Consequently, cq and Cp can be classified as 
fast functions as compared to the more slow functions no and Up. Then the system (18) of the equations for the fast 
functions can be approximately solved by keeping the slow functions, no and Up, as quasi-integrals of motion. From 
Eq.(18) we get the equations 



rf^co , ., . ,rfco 



df^ 



dt 



dt^ 



.,„ + A.,f 



— anp{anQ ~ Alo) 



— ano{anp + Aw) 



co = 0. 



(22) 



with the initial conditions (20) and 



co(0) 



cp{0) 



-(3* 



(23) 



where the dot means time derivative. The solution of Eqs.(22), with no and Up kept fixed, writes 



Co 



fit Q!(no — n„) 

r 2 +^ ^ 


- Aw 


sin- 




exp j 


-\^-- 



exp<, --{a- Au))t\ , 



(24) 



with the effective Rabi frequency given by the expression 

Q, = [Q;(no — rip) — Aw]^ 



(25) 



Then for the slow functions wc obtain 



|/3|2 . 2 m 

no = l-^sm -, 



f}2 "'" 2 • 



(26) 



The functions in Eqs.(26) describe the time evolution for the population of the ground state and of the chosen excited 
state. The form of these functions is similar to that one meets considering the Rabi oscillations [17]. However, it 
is worth emphasizing that, contrary to the linear case which can be recovered by putting a = 0, the expressions in 
Eqs.(26) are, actually, the equations for no and Up since the effective Rabi frequency (25) itself depends on these 
populations. Because of this, the solution of the equations from (26) will not result in simple sinusoidal oscillations. 
Consider, for example, the case when the detuning is such that it satisfies the relation 



Aw = 0. 



(27) 



Then Eq.(25) gives 



VL 



Aa^r 



I/3P- 



(28) 



In that case Eq.(26) shows that the ground-state level becomes empty while the upper resonant level completely 
populated, i.e.. 



no(tfc) = 0, np{tk) = 1, 



at the moments of time 



*fc = ^(l + 2fc) (fc = 0,1,2,...). 



(29) 



As far as no -^ 0, when t -^ tk, then the effective Rabi frequency (28) softens, il -^ \f3\, and the motion around 
t ^ tk slows down. Hence, the system spends more time on the upper level than in the ground state. And if at the 



moment t — tk we switch off the pum.ping field (7), then we shall get an inverted system with all atoms being in the 
non-ground state. Another way of obtaining an inverted system could be by adiabatically varying the detuning, as 
in the regime of adiabatic passage [17], till we reach the compensation condition (27). The latter, certainly, has sense 
only when the detuning continues to obey the resonance conditions (10)- (12). 

If the compensation condition (27) cannot be satisfied, then it is impossible to transfer all atoms from the ground 
state to the chosen excited state. However, it is always possible to populate these states equally. Really, consider the 
resonance case, when Alo ~ 0. Then at the moments 

i;;-2^(l + 8/c) (fc- 0,1,2,...) (30) 

we have 



1 

2' 



that is, both states are equally populated. 

It is worth paying attention to the following. The characteristic frequency of solutions (24) and (26) is the Rabi 
frequency (25). Looking back at Eq.(14), we see that in general, there should exist as well solutions with extra 
frequencies being combinations of basic frequencies entering Eq.(14). One may ask a question "when such extra 
frequencies could appear?" 

Recall that solutions (24) and (26) correspond to the guiding centers which constitute the main approximation in 
the method of averaging [21] and in the guiding-center approach [22]. The general form of the population amplitudes 
c„ is defined by Eq.(14). Denoting the general solution to the latter equation by cff", we may present it as a sum 

in which c„ is given by the guiding centers in (24) and (T„ is an additional ripple oscillating around the guiding centers. 
The characteristic frequency of the latter is the Rabi frequency (25). The ripple solution ct„ can, in turn, be written 
as a sum of terms with characteristic frequencies that are essentially higher than the Rabi frequency. In this way, the 
guiding center c„ represents the main harmonic while the ripple solution (T„ represents a sum of higher harmonics. 
Averaging the general solution c^*^" over the largest characteristic period corresponding to the higher harmonics gives 
the guiding center c„. Therefore, being interested in the average behaviour of solutions, one accepts the guiding center 
as the main approximation. Moreover, the ripple term (t„ not only oscillates much faster than the guiding center but 
the amplitude of the former is smaller than that of the latter. 

In order to concretize what is said above, let us substitute the general solution c-^^'^" = Cn+<Jn into Eq.(14). Introduce 
the notation 

A V 

-^^mjn ri * nin 

0^7njn ^ T 7 Pmn ^ T ■ 

Using equations in (18) for the guiding centers, we obtain the equation 

for the ripple term. Here c„ are the guiding centers defined by Eqs.(18) and (24). As is evident, the equation for the 
ripple term contains various higher harmonics, as a result of which the ripple solution a^ oscillates faster than the 
guiding center. 

Now, let us explain why the amplitude of the ripple term is smaller than that of the guiding center. Introducing 
the notation e = max{a,/3}, we see from (25) that the Rabi frequency $7 ~ e. Taking this into account and looking 
at (24), we conclude that the amplitude of the guiding center c„ is of order unity, c„ ~ e/fi ^ 1. 

The ripple solution (t„ can be presented as a sum of harmonics with amplitudes of order e/fJ^, where il^ is a 
characteristic frequency of a z^-harmonic. Among these characteristic frequencies there are various conbinations of uj 
and LUmn^^- Remember that, according to the quasiresonance condition (10), we have uj ~ Wpo ~ ^- Consequently, the 
characteristic frequencies of harmonics are of order O^ ^ ve, with i^ = 2, 3, . . .. Thus, the corresponding amplitudes 
are of order s/fli, ~ l/v. Hence, the amplitudes of higher harmonics diminish as l/iy with increasing i^ = 2, 3, . . .. 

In this way, we see that the guiding centers in (24) really constitute the main approximation to Eq.(14). In this 
approximation, the ripple solution, which oscillates faster and has smaller amplitude, can be neglected. If needed, the 




higher-order harmonics can be taken into account by means of perturbation theory. Such situation is common for 
the method of averaing and guiding-center approach [21,22]. 
The space-time distribution of atoms is given by the density 

P(r,t) = Poi^,'t) + Pp(r,t), 
in which 

p,(7,t)=iVn,(i)|^(7)|2 
is a partial density for j ~ 0,p. These densities, are normahzed to the total number of atoms 



N ^ I p{r,t)d r 
and, respectively, to the number of atoms 

TV,- = / pj{r,t)d r= Nnj{t) 



in the corresponding states. Since the spatial dependence of the wave functions for different states is different, we 
may get condensates with different space distributions. In general, such condensates will coexist, though, if the 
compensation condition is achieved, a pure non-ground-state condensate can be realized. 

In our consideration we have assumed that the system is initially cooled down so that all atoms are condensed in the 
ground state. The possible admixture of non-condensed thermally excited atoms has been neglected. Such a picture, 
as is known, is admissible for sufficiently low temperatures below the condensation point. If the temperature is kept 
low during the process of the resonant pumping, we may continue disregarding thermal excitations. Their role becomes 
important only after we switch off the pumping field. Since, during this pumping, the system has acquired additional 
energy, the latter can be redistributed among the energy levels through some relaxation mechanism. The role of such 
mechanism will be played by the interactions between condensed and non-condensed atoms. The thermally excited 
atoms will form a kind of a heat bath providing the possibility of relaxation to equilibrium. 

III. STATIONARY STATES 

Stationary states for the nonlinear Hamiltonian (3) are defined by the eigenvalue problem (6). The confining 
potential, typical of magnetic traps, is well described by the anisotropic harmonic potential 



m , 

y 



U.i7)^^iulx'+c.y+ulz'). (31) 



It is convenient to pass to dimensionless quantities measured in units of the characteristic oscillator frequency loq 
and length Zq given by the expressions 



LOO = (uJxiOyUJz)^^^, lo = \ 

V muo 



mujo 
The anisotropy of potential (31) is characterized by the anisotropy parameters 



(32) 



Ai = — , A2 = ^, A3 = — . (33) 

Wo 1^0 "^o 



Define the dimensionless coordinates 



x\ = -r, X2 = — , X3 = — , (34) 

to to '0 



forming the vector x— {xi,X2,X3}. The dimensionless interaction parameter is 

mA , , a 



g^-2-==47r(7V-l)-. (35) 

fi lo h 



Introducing the dimcnsionlcss Hamiltonian and wave function, respectively. 

g(y 

htoo 



H^^, ^(x) ^ /o^/^r ), (36) 



we obtain 

H-lt{-^ + ''-')+9m'- (37) 

Even when the scattering length a is much less than the oscillator length Iq, the interaction parameter (35) can be 
very large because of the great number of particles N. This situation is similar to that existing for large clusters [23]. 
Actually, a group of atoms trapped in a confining potential also forms a kind of a cluster. With a large interaction 
parameter, one cannot apply the standard perturbation theory in powers of g for calculating the eigenvalues of the 
Hamiltonian (37). Nevertheless, one may employ the renormalized perturbation theory [24-30] which, as has been 
shown by a number of examples [14-16, 24-31] successfully works for arbitrary values of the coupling parameter, as 
well as for all energy levels, providing a good accuracy with the maximal error around 1%. The first step of this 
approach is to choose an initial approximation containing trial parameters that, in the following steps, will be turned 
into control functions controlling the convergence of the procedure [24] . These control functions are to be defined 
from the fixed-point and stability conditions [29,30]. One of the simplest forms of the fixed-point condition is the 
minimal-sensitivity condition [28] which, for the first-order approximation, is equivalent to the variational condition 
for an energy functional [27-32]. 

As an initial approximation, we may take the Hamiltonian 

i—l ^ * ^ 

in which the effective frequencies Ui, with z = 1, 2, 3, will be control functions. The eigenfunctions of the Hamiltonian 
(38) are the oscillator wave functions 



V'«(a;) = JJV'«,(a;j), n = {ni,n2,n3}. 



i=l 

where n^ = 0, 1, 2, . . . Respectively, the eigenvalues of Hq are given by 

3 



4"^=E"^("^ + 0- (39) 



Perturbation theory is to be constructed with respect to the perturbation AH = H — Hq, which is 

3 

2 



The eigenvalues, in the first-order approximation, are 

i?W(A,.9,«) = i?(")+Ai?„, (40) 

where, for compactness, the notations A = {Ai, A2, A3} and u = {ui, U2, U3} are accepted and 

A£;„ = (V'„,Ai7Vn)- 

The control functions ui = Ui{X,g,n) can be found from the variational condition 

^ii;«(A,5,w) = 0, (41) 

which is a simple form of the fixed point condition. Substituting the found Ui into Eq.(40), we have 



For what follows, it is useful to introduce the notation 



(42) 



9n — 9^n-) 'Jn — I I '-^nit 



i^l 



in which 



^rii 



d^nJM^nJ') 



y^I^R 7r(2".n,!)2 



+00 



i/^^(x)exp(-2x2)dx, 



where Hn-{x) is a Hermite polynomial and Ui{n) = Ui{X,g,n). 
Eq.(40), with notation (43), can be written as 



El^H^g.u) == -^ ('^»+ 9 ) ("» + — ) + \/'"lW2W3 



(43) 



(44) 



Condition (41) results in the equation 



"-i + 2 ) ("i ^ ^i ) + U2y/UiU2U3 gn = 0. 



Using Eqs.(44) and (45), for the spectrum in Eq.(42) we have 



en(A, 5) = - ^ f ni + - j f u, + 5 



i=l 



\2 



(45) 



(46) 



where m — Ui{X,g,n) arc defined by Eq.(45). 

To understand better the properties of the spectrum (46), let us consider the weak and strong coupling limits. In 
the weak-coupling limit, when gn — > 0, the solution to Eq.(45) writes 



VA1A2A3 , / , Ej=i ('^j + 5) Aj - (ri» + i) A, ^ A1A2A 
Ui — Ai — -n \\9n + \ -^i 



2 (n, + i) 
Substituting this into Eq.(46), we get 



mUi^^ + h) 



4(r., + i)^A 



3 2 

9n- 



(A, g) =i X! ( "•» + 9 ) ^* + V-^i^aAs 5n - o X! 



A1A2A3 



2 



as 5„ ^ 0. 

In the strong-coupling limit, when g„ -^ 00, the solution to Eq.(45) reads 

„. ^ ("^ + l) ^' „-2/5 

1/5 ^" 



nL(n. + i)A| 



And for the spectrum (46) we obtain 



e«(A,5) ~ - 



nf".4)-^: 



i=i 



2 

2 / - J 



1/5 



i2/5 



(47) 



(48) 



(49) 



(50) 



as gn -^ 00. 



IV. GROUND STATE 



The ground state plays a special role for the phenomenon of Bose condensation. Therefore, we pay a little more of 
attention to the case of n^ = 0. Since Jq = {2n)^^^^, the effective interaction strength (43) becomes 



50 



9 



(27r)3/2- 
Eq. (46), defining control functions, simplifies to 

u'^ + 2goUiy/uiU2U3 - \f = 
And the spectrum (46) reduces to 

3 



eo(A,5) = ^E("^ + 5|)- 



(51) 



(52) 



(53) 



In the weak coupling limit, when go -^ 0, Eq.(52) yields 



Ui — Ai - VA1A2A3 go + 



Ai / Y^ \ \ I ' A1A2A3 



vJ=i 



A,; 



9l 



And for the spectrum (53), we have 



eo(A,5) ~ -(Ai + A2 + A3) + VA1A2A3 50 - -j-^(AiA2 + A2A3 + A3Ai)5o, 

as 50 -^ 0. 

In the strong-coupling limit, when go ^ 00, for the control functions we get 



^» 0-2/5 

(2AiA2A3)2/5^0 ' 



and, respectively, for the spectrum we find 



eo(A,g)c.-(2AlA2A3)2/^5o^^ 



(54) 



(55) 



(56) 



(57) 



For an axially symmetric potential, with Ai — A2 7^ A3, the strong coupling limit (57) reduces to that found by Baym 
and Pethick [33]. 

For an isotropic potential, for which Aj = 1, there is only one control function given by the equation 



Then the ground state energy (53) becomes 



+ 250"^/^ -1 = 0. 



4 V u 



In the weak-coupling limit, when go ^> 0, Eq.(58) gives 



And the spectrum (59) is 



u~ l-go + 2go- 



Mg) -2+90- -^90, 



(58) 



(59) 



(60) 



(61) 



as 50 -^ 0. 

In the strong-coupling limit, as 50 ^ 00, we have 



10 



u 0, (2.go)-'/^ (62) 

and, respectively, 

eo(5)^^(25o)'/^ (63) 

For the atoms with negative scattering lengths, as in the case of ^Li or ^^Rb, the interaction parameter (35) is 
negative. All weak-coupling expansions, such as Eqs.(48), (55) and (61), are the same for the negative g -^ —0. When 
g increases, the real solution for the spectrum exists till some critical value gc after which it becomes complex. Thus, 
for the isotropic potential, the stable ground state is defined for gg < 50 < 0, with 

.9S - -^ = -0.267496. (64) 

From here, the critical interaction parameter is 

gc = (27r)3/2gg = -4.212960. (65) 

This, according to notation (35), defines the critical number of atoms 

N^^l + ^ (66) 

that can condense in a stable state. After g < g^, the spectrum becomes complex. The latter means that the 
corresponding states are not any more stationary but quasistationary. The lifetime of a quasistationary state is 
defined [34,35] as 

'Tnig) = 7, — TT TTT- ^^^^ 

2ujQ\Im en{g)\ 

The complex spectrum of quasistationary states is defined, as earlier, by Eq. (46), with Ui from Eq.(45). These 
equations, in the case of negative (/ < 0, have several complex solutions, of which we must choose that continuously 
branching, at g — gc, from a solution that is real at g ^ —0. Solving these equations for the isotropic ground state, 
we find 

i?e u ~ 0.705470g"2/^ + 3.850396g"^/^ + 12.402511.g"^ 

Imu~ 2.171212.9^2/5 _^ 2.797476.g"*^/^ (68) 

for g ^ 1. For the real and imaginary parts of the energy we obtain 

Re eo{g) ~ O.WndSg^^'^ + 0.529102.g"2/5 _^ 1.4438999-6/5 + 31.0062775^^ 

Im ea{g) ~ -0.520739^2/5 + 1.628409g"2/5 ^ i.049054g-6/^ (69) 

Therefore, for the lifetime (67) we have 

, , 0.960174 
Mg) ^ 57^. 70 

If the number of atoms, with a negative scattering length, exceeds the critical number given by Eq.(66), then only Nc 
atoms can form a stable Bose condensate, excessive atoms being expelled out of the condensate during the time (70). 

V. EXCITED STATES 

The energy of any excited state is defined by Eq.(46). In order to illustrate the relative to each other behaviour of 
these states, let us write explicitly the corresponding formulas for several first levels. Consider, for the isotropic case, 
the energy levels 6100(5), ^110(5) and 6200(5)- For these levels, the integrals J„. entering notation (43) are 

J - ^ J - ^ J - ^^ 



V27f' 4V27f' 64^/2^' 
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In the weak-coupling limit, the energies of the first three excited states, compared to that of the ground state, 
behave as 

eo(g)~^ + 0.063494g, 

eioo(5)^^ +0.0476205, 

eiio(ff)~^ +0.0357155, 
6200(3)^^+0.0406765, (71) 



when 5^0. And in the strong-coupling limit, we get 

eo(5) ^ 0.547538ff2/5, 

eioo(5) - 0.607943<?2/5, 

e2oo(g) ^ 0.632193.92/5, 
eiio(5) - 0.6750125^/5, (72) 

as 5 ^ 00. 

Notice an interesting and important fact that the energy levels 6110(5) ^^id 6200(5) cross each other as 5 varies. 
Really, from Eq.(71) it follows that 6110(5) < ^200(5) at small 5, while Eq.(72) shows that 6110(5) > e2oo(,9) at large 
5. This crossing of levels demonstrates that the latter cannot be classified being based on the harmonic oscillator 
spectrum. To correctly classify the levels, one has to calculate their energies for each given interaction strength 5. 
At the same time, the strong distortion of the harmonic-oscillator spectrum plays a positive role for our purpose, 
making the spectrum of interacting atoms nonequidistant, which is necessary for the possibility of transferring the 
ground-state Bose condensate to a non-ground state, as is explained in Sec. II. 

The transition frequency, measured in units of luq, is given by the difference 

(^n{g)=en{g)-eo{g). (73) 

For the low-lying excited states this leads to 

ujioo ^ 1 - 0.0158745, 

WHO ^ 2 - 0.0277795, 
W200 ^ 2 - O.O228I85, (74) 

in the weak-coupling limit, 5 ^ 0, and to 

wioo(5) ^ 0.0604055^/5, 

c^2oo(5) - 0.0846555^/5, 

wiio(5) ^ 0.12747452/5, (75) 

in the strong-coupling limit, as 5 -^ 00. Eqs.(71)-(75) show that sufficiently strong interaction parameter 5 makes 
the spectrum nonequidistant, which favours the possibility of the resonance pumping of atoms from the ground state 
to an excited level. 

The shape of the atomic cloud, in a given state, can be characterized by the aspect ratio 

12 



R 



13 



{A)^"' 



(76) 



in which {x^) is a mean-square deviation in the x-direction, and (X3) is a mean-square deviation in the z-direction. 
Similarly, we may define the aspect ratio i?23- In the case of cylindrical symmetry, i?i3 = i?23- 
The aspect ratio (76) can be written as 



i?i 



ni 



W3 



(na + i) Ml 



1/2 



(77) 



In the weak-coupling limit this yields 



''"-V(-3 + i)Arr + I(^^0^ 



1 



1 



(ni + i) Ai (n3 + i) A3 
when g ^f Q. For example, for the ground state, Eq.(78) reduces to 

^3 



n 



13 



1 



VA1A2A3 /I 1 



2(27r)3/2 VAi A3 



In the strong-coupling limit, Eq.(77) gives 



i?i3 ~ — {g^ 00) 

Ai 



(78) 



(79) 



(80) 



for all energy levels. 

For the ground state, our results agree with the aspect ratio found by other authors [33,36,37]. And for excited 
states, we predict an interesting fact that the aspect ratio in the strong-coupling limit becomes asymptotically, as 
g ^ 00, independent of level numbers, according to eq.(80). 

VI. INTERACTION AMPLITUDE 

An important quantity defining whether it is possible to transfer all atoms from the ground state to an upper level, 
or only half of them, is the interaction amplitude (16). For convenience, let us define the dimensionless amplitude 



a-nig) 



Wo 



(81) 



which may be presented as 



where Ui{n) = Ui{X,g,n) and 



a„(g) = a/ui(0)m2(0)u3(0) /„c 



/n = JJ^,, n= {ni,n2,n3}, 






0^ ^^2' 

The last integral can be expressed as 



"'"^j! J-00 I Ui{n) J 



_ (-l)"»r(n, + l) l 2Q-l ^f _ l__ ^ 
"•" ^nJCr V 2C. ^ ""'' "''2 ''"^' 
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through the gamma function r(. . .) and the hypergeometric function F{. . .), with 

Ui{0) + Ui{n) 



Q 



In the weak-couphng hmit, using the equahty 



2w.(0) 



F ( -rii, -n,; - - n^; 1 ) = (-1)"' 



we find 



an{g) 



V2A1A2A3 

47^3 



'j^ r(n, + i) 



(82) 



3/10 



as 5 -^ 0. And in the strong-couphng hmit, we obtain 

a„(g) 0, (AiA^Ag)^/^ f ^j x 

X n ^"'^" 






where g ^ 00 and 



- r(n, + iU^^f (^-n„-n,;i-n,;C 



^'=2 + 



i2/5 



(83) 



27r 



!nli(^. + i)^^. 



1/5- 



To compare the interaction amphtudes for the first several levels with the corresponding energies and transition 
frequencies studied in Sec.V, let us take the isotropic case. Then, in the weak coupling limit, Eq.(82) yields 

aooo(5)-0.06349g. 



"100(5) - 0.03175g, 

aiio{g) ~0.01587g, 
^200 (3) - 0.02381.g, 
as g ^ 0. And in the strong-coupling limit, from Eq.(83) we derive 

aooo(3)- 0.2190252/5, 

2/5 



(84) 



aloo(ff)^0.183045"/^ 

anoi9)-0.U759g^/^ 
a2oo{9) - 0.1755952/5, 



(85) 



as 5 — > 00. 

Comparing the interaction amplitudes (84) and (85) with the transition frequencies (74) and (75), we see that the 
compensation condition (27) could be accomplished in the weak-coupling limit or in the intermediate region of g, 
where an{g) <C w„(g). In the strong-coupling limit, an{g) becomes of the order of uJn{g) and condition (27) cannot 
be satisfied, since \Au>\ <C LUn{g)- This means that, with increasing g, that is, with increasing the number of atoms in 
the ground-state condensate, it becomes more difficult to transfer all these atoms to an upper level. For large atomic 
clouds, only a half of the condensed atoms can be pumped up to an upper level. 
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VII. DISCUSSION 

We suggested a mechanism for creating nonequilibrium Bosc condensates in non-ground states. The possibihty of 
transferring either all atoms or only a part of them to a non-ground state depends on the parameters of the system. 
In principle, these parameters are changeable, and we think that it is possible to adjust them so that to effectively 
realize the mechanism suggested. For concreteness, let us look at the parameters typical of magnetic traps used now 
for condensing atoms in the ground state. 

In the JILA trap [1,38] the atoms of ®^i?6, with the mass m ~ 1. 445 x 10^^^ g and the scattering length a = 
6 X lO^^cTTi, were condensed. The characteristic frequencies of the confining potential are oj^ — ujy — 2tt x 132 Hz and 
Wz = 27r X 373 Hz, their geometric mean, defined in Eq.(32), being ujq — 2Tr x 187 Hz. The anisotropy parameters 
(33) are Ai = A2 = 4r and A3 = 2. The oscillator length Iq = 0.788 x lO^^'cm. Therefore, the interaction parameter 
(35) is 

5 == 9.565(iV - 1) X 10^2^ 

For the number of atoms A^ ^ 10'^, one has g ^ 100, which corresponds to the strong-coupling limit. 

In the MIT trap [3,39] a condensate of ^^Na, with m ~ 0.382 x lO^^^g and a ~ 5 x lO^^cm, was realized. The 
characteristic frequencies are Ux = ujy = 'in x 320 Hz and w^ = 27r x 18 Hz, so that their geometric mean is 
Wq = 27r X 123 Hz. The anisotropy parameters are Ai = A2 = 2.610 and A3 — 0.147. The oscillator length is 
lo — 1.890 X IQ^^cm. Thus, the interaction parameter (35) is 

g = 3.324(A- 1) X 10"^. 

For iV ^ 5 X 10^ — 5 X 10®, this gives g ^ 10'' — 10^, which certainly is in the strong coupling limit. 

In the RQI trap [2,40] the atoms of ^Lz, having m = 0.116 x 10~^^g and the negative scattering length a = 
— 1.5 X lO^^CTO, were cooled down to the Bose condensation regime. With the characteristic frequencies lUx — 
2tt X 150.6 Hz, ojy = 2tt X 152.6 Hz, and w^ = 27r x 131.5 Hz, the geometric-mean frequency is utQ — 2tt x 145 Hz. 
The anisotropy parameters (33) are Ai = 1.039, A2 = 1.052, and A3 = 0.907. The oscillator length Iq = 3.160xlO^'*c?7i. 
The interaction parameter (35) becomes 

g = -0.597(iV-l) X 10^2. 

Since the confining potential is almost isotropic (Ai w A2 ~ A3) we may take the value (65) for the critical interaction 
strength allowing yet a stable condensate. Then, the critical number (66) of atoms allowed in the condensate is 
Nc = 707, which is in agreement with the estimates of other authors [33,37] giving Nc ^ 10'^. 

In this way, the trapped clouds of *^i?& and '^^Na corresponds to the strong-coupling limit, and that of ^Li, to 
an intermediate regime, with g < 4. Therefore, if the number of condensed atoms in a cloud is iV > 10'^, it would 
be difficult to transfer all of them to an upper level. However, as is discussed in Sec. II, it is always possible to 
transfer a half of the ground-state atoms to a non-ground-state level, thus, creating two coexisting condensates of 
the same species in two different quantum states. In addition, it is, probably, possible to construct specially designed 
traps allowing to pump up all atoms from the ground state to another level, forming a pure non-ground-state Bose 
condensate. In any case, both these possibilities, being realized, would produce a new kind of quantum ultra-cold 
matter which may reveal interesting and unexpected properties. 
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